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TECHNICAL MEMORANDUM X- 64964 


ELASTIC DEFORMATION OF THE GYRO ROTOR 
IN THE RELATIVITY GYRO EXPERIMENT 

INTRODUCTION 


The elastic deformation of a spinning sphere is an old problem which was 
first solved by Chree in 1889 (IJ. This solution is needed to evaluate the elec- 
trical torques on the rotor in the Stanford Gyro experiment. In this experiment, 
the elastic deformation due to spin will produce nonradial forces on the ball 
which will lead to drift-producing torques. These torques must be understood 
and calculated to estimate the ultimate accuracy of the experiment. 

In th" original derivation by Chree [11 and in subsequent treatments [2, 

3] , the results are exhibited without a derivation. A simplified derivation which 
points out the assumptions which are made in the derivation is presented here. 
Spherical coordinates are used so that the end result can readily be used as an 
input to the electrostatic torque calculation. Numerical results are presented 
for the most likely configuration in the gyro experiment, a quartz ball spinning 
at 200 Hz. The fractional change in difference of principal moments of inertia 
A 1/ 1 is then calculated including density changes in the ball. This result is 
needed for the gravity gradient torque calculations. 


■ DEFORMATION CALCULATION 


The problem is attacked by solving the Navier equations in elasticity 
theory with the appropriate boundary conditions. These equations must be 
written and solved in spherical coordinates because for our electrostatic 
torque calculations, the shape of the ball must be developed as a Fourier 
cosine series in the polar angle 0 . To verify that the shape contains only an 
a .2 cos 20 term and to calculate a 2 , we must work in polar coordinates. We 
start from first principles to bring out any approximations or assumptions 
made in the derivation. This presents a problem since most books give the 



equations only in cartesian coordinates. The simplest way to find the equations 
in spherical coordinates is to write the equations in their general form with 
covariant derivatives present and then find the form of the covariant derivatives 


in spherical coordinates. To begin» one calculates the strain tensor 
terms r C the displacements using 


in 



where D/ Dx 


i 


denotes covariant derivatives. 


When this is done and the identi- 


fication of the physical components of u^ 
equation 48.17 of Sokolnikoff [4]. 



is made, one can obtain 


In spherical coordinates, the line element and Christoffel symbols are 


ds^ = dr^ + r^dO* + sin^ dd0^1 
r 33 = -sin 0 cos 6 Tii = -r rjj = -r sin* $ 
rit = rjg = rli = cot e . 

Identifying the ’’physical components” j 5) of u. as u =uj, u =U 2 /r, and 

i r Q 

u = U 3 / r sin 0 , one finds for the physical components of e.. from equation ( 1 ): 
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( cont' d) 


where the physical components of e.. are J 's/g^e . The dilatation 
is given by ^ 
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The equations of equilibrium in tensor notation are given by 


.. Dt.. 

^ F. = 0 

Dx 


(3) 


where t. is the stress tensor and F. is the body force vector, in this case 
ij 1 

representing centrifugal force. The first term in this equation must be cal- 
culated in spherical coordinates. Using the definition of covariant derivative 
and the Christoffel symbols, the following equations are obtained: 
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( 6 ) 
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« -F 


These equations contain the physical components of r . and agree with equation 
48 . 18 of Sokolnikoff [ 4] . 

To obtain equations containing only the u^* s (the Navier equations) , 
the s must be expressed in terms of the u^' s . This is done by the 
generalized Hooke’s Law valid for a homogeneous Isotropic body 


Ty = xagy * a.e,j 


where X and n are the Lam4 constants. Using equation (2) , we obtain for 
the physical components of : 
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= xe + ^ 


where we have assumed that the u * s are independent of (and u o) . 

i 

These equations can now be 'ubstituted into equation ( 4 ) , from v4iich we obtain 
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Substituting equation (7) into equation (5), we obtain 


r 90 r 


Q 2 

“T (ru ) - ^ ^ = -F 

8r* ' d' r dr$e 0 


( 9 ) 


Equations (8) and (9) are the Navier equations of elastic equilibrium written in 
spherical coordinates with no <f> dependence. They agree with equation 110* 
of Honeywell’ s report No. 20831 FR [ 6J . 

They must be solved for the case of a body force given by the centrifugal 
force on the spinning ball. The solution must satisfy the boundary condition 
that the surface of the spinning sphere being considered is free from stress. 
That is, at r = a , we must have 


T 

rr 


0 , 



The body force can be (terived from a scalar potential given by 


( 10 ) 


V 



w*r* -^in* 0 
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where p is the density of the ball and w its angular velocity. Then, in 
spherical coordinates. 


F = “ — = pcu*r sin 0 
r or 

( 11 ) 

1 av 2 

" TT ~ Pw r sin 0 cos 0 


Now the solutions to equations (8) and (9) are written with the body force given 
by equation (ll). This solution sattsfies the boundary conditions of equation 
( 10 ): 
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It can be verified that this is a solution of equation (8) as follows. First, if 6 
is calculated from this solution, 


(X 2p) ^ = pw^a* |(x + 2p) |^P2(cos 0)(5ai + 2/52) 


2r r 
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The otlu3r terms in equation (8) give 



Adding U^se three terms and using the identities 

q!i(5X+7^) + Pi{2\) - ^ 

10(X + 2^) aj 

V 


we find that the left side of equation (8) gives -pui^r sin^ 6 in agreement with 
equation (11) . The three terms in equation (9) are 



Combining these terms, we verify that the left side of equation (9) gives 
-pcu^r* sin 0 cos 0 in agreement with equation (11) . Therefore, this verifies 
that our solution satisfies equations (8) and (9) . It remains to show that the 
solution also satisfies the boundary conditions in equation ( 10) . Usiiq; 
equation (7) with the differentiated form of equation (12) , we find that at 
r = a 


7 



^rv “ P‘*>*®*{P2(coB «)lof|(6X+ 6fi) + 2X^l + 2^a|] 
+ laj(5X+6fi) + a4(3X+ 2^)]} 


The second square bracket is obvi<^ly zero and, using the expression for , 
fii and 0 ^ 2 , t)% first square bracket is also zero. Again, using equation (7) 
and t)« differential form of equation (12) , we find that at r « a 


^rO ~ ® 0(^3(xi - 3of2 + 2^i) 


and this is also zero by the definition of a^, atf and . This (temonstrates 
that equation (12) satisfies both botmdary conditions in equaticm (lo) . 

Evaluating the solution at 
2oojV ( _ , 

“r ° 

where we have used the definitions of the Lam4 constants in terms of Young's 
modulus and Poisson* s ratio: 

X 

^ 2(X + ^) 
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From this we note that the fractional difference A 
rial and polar radii is 


between the equato- 


A 


= f) - u^(9 . 0)] 
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We also note that since 


3 /l + cos 2fl\ 1 
P«(co.«» - J j-- . 

the coefficient of the cos 20 term, which we will call a 2 , is given by 



Thn ccH^fficient &2 i^cessary to calculate electric torques on an Clectro- 
s^.atic i^^ro since they ' .re calculated in terms of harmonic coefficients of 
which a; is the lovest. 

Now we can compute A and &2 numerically for a quartz ball spinning 
at 200 Hz. We will use the following parameters: 

p = 2,2 g/ cm^ 

0 } = 200 cps 

a = 0.75 in. » 1.9 cm 

E = 10.4 X VD* Ib/in,* « 7. 17 x lo“ dynes/cm* 

t' * 0.16 

= 1.76X 10"‘ 

A = 5. 62 X 10“* 

aj = -2. Ip in. =* -5. 33 x lo“* cm 
Both these quantities scale as . 
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MOMENT OF INERTIA CALCUIATION 


Using the results of the last section, we can calculate the difference in 
the principal moments of inertia produced by the centrifugal expansion. This 
calculation is needed to evaluate the effects of gravity gradient torques on the 
gyro experiment. In the expression for gravity gradient torques, the qusuitity 
A 1/ 1 is needed, where 


AI = Ij3 - fdu^In) = I / dvlr^ - 3z*lp(r) (14) 


and I|i, I22, and I33 are the principal moments of inertia. There are two con- 
tributions to Al/ 1: the first due to the oblateness of the ball, and the second 
due to density changes in the ball. The first effect is calculated as follows: the 
expansion of the ball given by u {&) produces an extra infinitesimal element 
dm given by 


(im = pu^( 0)(2 7ra sin 0) (ad 0 ) 


and one then needs to calculate 


AI 



a^[l - 3 cos^ 0 j dm 



P2(cos 0)dm 


where P2 is the second Legendre polynomial. Using equation ( 13 ), this 
becomes 
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The second integral is zero because of the orthogonality of Legendre polynomials 
( Pq - 1) , and the first int^n^al is equal to 4)r/5 by the normalization int^prals 
for P 2 . Hence, 


AT _ pj 

15E (7 + 5p) 


and, since 1= 8/15 Trpa®, we finally have 


A1 pcj^a^ ( 1 + v) {2+ v) ( 5X + 4u \ 

I " E (7+ 5v) * 2p \19X+14pj 

which is equal to A . Therefore, we have the result that for the distorted sphere 
of uniform density, AI/ 1 = A [ 3] . However, in actuality, there will be a change 
in density when the gyro is spun up which will also produce a contribution to 
AI/ 1 . This is given by 

7T 2 

AI' = -2ir f sin 0 d0 J r^ dr ” U - 3 cos d]pQ , 

0 


where the integral is now over the entire volume of the rotor ( assumed spherical) 
and -0 is the fractional density change Ap/p . Using the expression for 0 
computed from equation ( 12) , tliis integral becomes 


^ 0 


f [P 2 (cos 0 ) 1 * Sin 0 d0 
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’3 (7 + 5t;) a^ 


+ f lp 2 (cos 0)) sin 0 do 
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i Ll. t J'l 

3 (1 - 
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As before, the second 6 integpral is zero and the first is 2/5; therefore, this 
expression gives 




-47rp^coV ( 1+ t/) ( 1 - 2i/) 2 

21E (7 + 5t/) 5 


AP _ -0(J^a? ( 1 + t^) ( 1 - 2t^) 

P ^ 7E (7+5t/) 


-PcuV 1 
7 19X+14p 


This gives for the total change in AI/ 1 

AI ^ ogjV 35\ + 2 6u 
I 14 p 19\+14p 


which is agreement with Reference 3. This gives a reduction in AI/ I of 


1 - 2v 
7{2+ v) 


4.5% 


Therefore, for the quartz sphere at 200 Hz, we have 


Y = 5. 37 X 10"® 


While this is a small number, the usual expression for gravity gradient torques 
involving AI/I indicates that they will contribute one of the largest torques on 
the gyro in the final gyro experiment. For this reason it is important to know 
the actual magnitude of the AI/I term accurately. 
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